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CLAIMS 

WHAT IS CLAIMED IS: 

1. A method of computing FIR filter coefficients, comprising the steps of: 

inputting a filter order of a universal maximally flat FIR filter, a number of 
5 zeros at z=-l, and a parameter for a group delay at z=l, the filter order being a positive 
integer, the number of zeros being an integer equal to or more than zero, the parameter 
being a rational number; 

executing a first operation by a first recurrence formula which includes 
parameters for the filter order, the number of zeros, and the group delay, and provides 
10 coefficients in Bernstein form representation of a transfer function of the universal 
maximally flat FIR filter; 

executing a second operation by a second recurrence formula composed of 
additions, subtractions, and divisions by 2, by using a resultant of the first operation as 
an initial value; and 

15 extracting impulse response coefficients of the universal maximally flat FIR 

filter from a resultant of the second operation. 

2. The method according to claim 1, wherein: 

the first recurrence formula is expressed as 

b/ = (-l){(2d) bj.r + ( j - 1 ) bj. 2 '} / (N - j + 1 ) where 1 ^j^N with b 0 ' - 1 
20 andb-i , = 0, 

wherein the filter order is N, the parameter for the group delay is d, coefficients in 
Bernstein form representation of a transfer function of the universal maximally flat FIR 
filter are b/; 

the resultant of the first operation is expressed as B'={l,bi',...,b N -K\0,...,0}, 
25 wherein the number of zeros is K; 

the second recurrence formula is expressed as 

hi (p) = ( 1 + E ) hi*"" / 2 + (1 - E) hi_i (p " x) / 2 where 1 ^p^N, O^i^p with h 0 (0) 
= B 9 and h.^ = {0,.„,0}, 
wherein a sequence for computing impulse response coefficients of the universal 
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maximally flat FIR filter is expressed as hi to)s =(hy (p) ) s =(hi f o* ) ,hi > i (p) ,...), and an arbitrary 
sequence Ai is expressed as E* = E ( E j_1 A* ), E 1 Ai = EAi = Aj+i, E°Aj = Ai in which a 
forward shift operator satisfying the expression is E; and 

the impulse response coefficients extracted from the resultant of the second 
5 operation are expressed as hi = hi,o (N) where O^i^N 

3. A program for computing FIR filter coefficients, the program causing a computer to 
execute the steps of: 

determining every element of a single-dimension array B' using a filter order N 
being a positive integer of a universal maximally flat FIR filter, a number of zeros K at 

10 z=-l, K being an integer equal to or more than zero, and a parameter d for a group delay 
at z=l, d being a rational number, all of which are provided by inputs, by changing in 
sequence an index j from 1 to N-K in a recurrence formula B'fj] = (-l)X {(2d)B'[j-l] + 
G-l)B'D-2]} / (N - j + 1), the single-dimension array having N+l elements B'[j] where 0 
^j^N, in which an element B'[0] thereof is initialized to 1 and all the elements thereof 

15 except the element B ! [0] are initialized to zero; 

determining every element of a three-dimension array r by sequentially 
changing, in the order of indexes j, i, p, an index j from 0 to N-p, and an index i from 0 
to p, an index p from 1 to N in a recurrence formula r[p,ij] = ( r[p~l,i-l,j] - 
r[p-l,i-l J+l] ) / 2 + ( r[p-l,i,j] + r[p-l,i ) / 2, the three-dimension array r having N 3 

20 elements r[p,ij] where O^p^N, O^i^N, O^j^N, in which elements r[0,0,j] thereof 
where 0^j ^N-K are initialized to elements of the single-dimension array B'[j] where 0 
Si j N-K, and all the elements thereof except the elements r[0,0,j] are initialized to 
zero; and 

extracting elements r[N,i,0] of the three-dimension array r where O^i^N as 
25 the impulse response coefficients of the universal maximally flat FIR filter. 



